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ANALYSIS  OF  CONSTANT  FALSE  ALARM  RATE  SIDELOBE  CANCELLER  CRITERION 


L.  E.  Brennan  and  I.  S.  Reed 


I.V  INTRODUCTION 


'^In  this  final  report,  the  constant  false  alarm  rate  (CFAR)  detection  cri¬ 
terion  for  a  si delobe  canceller  (SLC)  system,  introduced  in  the  last  quarterly 
progress  report  is  found  completely  and  analyzed.  This  new  detection  test 
for  radar  exhibits  the  desirable  CFAR  property  that  its  probability  of  a  false 
alarm  (PFA)  is  functionally  independent  of  the  covariance  of  the  actual  noise 
field  encountered.  As  a  consequence,  such  a  CFAR  SLC  system  is  ideally  suited 
to  cope  with  the  newly  evolving  ^smart''^  jammer  threat  to  radar.  ^ 

W  I  — 

An  important  objective,  set  in  the  last  quarterly  progress  report^.X^» 
was  to  find  both  the  false  alarm  and  signal  detection  probabilities  of  this 
test.  The  first  and  most  important  of  these  two  goals  has  been  met.  The  prob¬ 
ability  of  a  false  alarm  (or  PFA)  of  this  CFAR  SLC  detection  criterion  is  de¬ 
rived  in  closed  form  in  this  report.  The  success  in  finding  the  PFA  is  due 

— 

primarily  tc  the  use  of  a  generalization  of  Cochran's  theorem^X^t^r'llT]. 


II.  FORMULATION  OF  THE  PROBLEM  / 

The  main  beam  of  the  radar  is  assumed  to  receive  a  complex  base-band  pro-  V 

cess  y(t),  which  is  sampled  to  yield  y(n)  for  (n  =  1,  2,  ...,  N).  Associated 

with  the  main  beam  of  the  radar  are  k  auxiliary  antennas  Xj(t),  called  AUXs,  -  — 

for  (j  =1,  2,  ...,  k).  Each  AUX  Xj(t)  is  assumed  to  have  a  gain  substan-  ^ 

tially  less  than  the  main  beam  gain— so  much  so,  that  it  can  be  assumed  that 

the  AUXs  contain  no  radar  echo  information  of  significance.  Each  AUX  x^  is  - 

J  ■  an/ 

sampled  in  synchronism  with  the  sampling  of  the  main  beam  in  such  a  manner  nt-s  c 


that  the  output  in  the  j-th  AUX  is  the  sampled  data  set  Xj(n)  for  (n  =  1,  2, 

....  N). 

To  restate,  let  |y(n)|n  =  1,  2,  ...»  n|  denote  a  set  of  quadrature  danodu- 
lated  and  sampled  main  beam  data  in  which  one  wants  to  locate  a  desired  radar 
signal.  Also  let 

jx(n)  =  [x^(n),  Xgin).  ...»  X|^(n)j^|n  =  1,  2,  n| 

denote  a  sampled  vector  set  of  the  k  AUXs  of  complex  base-band  data  which  is 
assumed  to  have  no  signals.  Next,  let  |s(n),  n  *  1,  2,  ...,  n|  be  the,  possi¬ 
bly,  phase  and  amplitude  coded  and  sampled  waveform  which  one  would  expect  to 

iS 

find  in  the  N  samples  of  the  main  beam  data.  Finally,  let  b  =  A^  be  a 
possible  complex  signal  envelope,  where  A  is  the  envelope  and  8  is  the  phase  of 
signal  s(n). 

The  radar  receiver  must  distinguish  between  two  hypotheses,  the  noise-only 
hypothesis,  denoted  by  Hq,  and  the  signal-plus-noise  hypothesis,  denoted  by 
If  Ngln)  denotes  the  noise-only,  sampled-data,  main-beam  process,  and  N(n)  = 
^Ni|nj . *'*k('’)]^  denotes  the  vector  of  AUX  data,  then 

y(n)  =  NQ(n) 

x(n)  *  N(n)  for  (n  =  1,  2,  ...,  N),  and 
y(n)  =  NQ(n)  +  b  s(n) 

x(n)  =  ^(n)  for  (n  =  1 ,  2,  . . . ,  N)  (1) 

are  the  outputs  of  the  main  beam  and  AUX  receivers  under  the  two  hypotheses. 

The  external  noise  field  is  primarily  due  to  time-varying  jamming,  not 
backscatter  from  clutter.  Thus,  it  can  be  assumed  that  the  noise  process 
[y(t),  x^(t),  ...,  X|^(t)]^ 


impinging  on  the  radar  and  its  AUXs  is  a  zero-mean. 


Gaussian,  vector  process  with  a  possibly  time-varying  covariance  matrix,  but 
Independent  from  time  sample  to  sample.  For  the  validity  of  the  latter  as¬ 
sumption,  time  sampling  needs  to  be  accomplished  at  the  Nyquist  rate.  The  CFAR 
SLC  detector  criterion  developed  In  the  next  section  Is  based  on  the  above, 
rather  general,  noise  model. 

III.  GENERALIZED  MAXIMUM  LIKELIHOOD  RATIO  (MLR)  TEST 

The  CFAR  sidelobe  canceller  test  of  the  following  section  uses  for  Its 
derivation  the  generalized  maximum  likelihood  ratio  (MLR)  test  originally  due  to 
Neyman  and  Pearson.  The  MLR  principle  Is  best  described  by  probabilities  on  a 
sample  space  X  with  a  parameter  set XI.  Let  x  £  X  be  a  point  In  sample  space  X. 
Also,  let  P(x;  d)  be  a  probability  density  function  of  x,  where  9  e  XI  Is  a 
parameter  vector  In  XI,  the  set  of  all  possible  values  of  6. 

Now  suppose  w  <  XI  Is  a  subset  and  XI-  w  Is  the  subset  of  XI  complementary 
to  w.  Next,  define  two  alternative  hypotheses,  Hq  and  Hi,  In  set  XI  In  the 
following  manner; 

«0 '  {®  “'  } 

H,  3  {9  €  a-  W  j .  (2) 

The  purpose  of  the  likelihood  ratio  test  Is  to  test  hypothesis  Hi  against  Hq, 
holding  the  probability  that  Hi  Is  accepted  when  Hg  Is  actually  true,  I.e., 
the  probability  of  a  false  alarm  (PFA)  below  some  fixed  error  level. 

An  Important  test  that  accomplishes  the  above  requirement  Is  the  maximum 


likelihood  ratio  (MLR)  test.  See  an  early  paper  by  Kelly,  Reed  and  Root  [3] 
for  more  details.  Let  K  be  some  fixed  threshold  and  define  the  likelihood  ratio 


Pfx*  fl) 

iM  , K  t  „ 

max  I 

0  €  w 

<  K  accept  Hq  ,  (3 

where  hypotheses  Hq  and  H-|  are  defined  In  Eq.  (2). 

If  &^(x)  are  the  maximum  likelihood  estimates  (MLEs)  of  9  for  a  given  x 
with  respect  to  hypotheses  |i  =0,  l),  then  L(x)  in  Eq.  (3)  is  equivalent  to 

P(x;®,(x)) 

“  P(x;  eo(x))  ‘  ^ 

In  this  form,  the  likelihood  ratio  L(x)  is  explicitly  related  to  the  W.Es  of 
the  parameter  vector  6  under  the  two  possible  hypotheses,  assuming,  of  course, 
that  such  estimates  exist.  The  MLR  test  is  used  next  to  find  an  adaptive  de¬ 
tection  algorithm  with  the  CFAR  property. 

Consider  the  (k  -t-  1)  component  complex  vector  of  main-beam  and  AUX  data, 

v(n)  *  I  “  [y(n).  •••»  \('^)]  * 

as  defined  in  Eq.  (1)  under  hypothesis  Hq.  Then,  under  hypothesis  Hq,  the  vec¬ 
tor  ^(n)  in  Eq.  (5)  has  the  joint  probability  density  function 


Po(i("))  ■ '(-‘"to  >  ■ 


(n)  M"^  ^(n) 


for  (n  =  1,  2,  ...,  N),  where  M  *  Ejv(n)  v*(n)|HQj  is  the  unknown  covariance 
matrix  of  v(n)  and  ||m  Ills  its  determinant.  Here,  the  asterisk,  denotes 


conjugate  matrix  transpose. 

The  covariance  matrix  M  in  Eq.  (6)  was  allowed  to  be  "slowly  time-varying" 
as  a  function  of  sampling  time  n.  Here,  M  is  allowed  to  vary  in  time  only  with 


-5- 


a  time  constant  which  somewhat  exceeds  N,  the  observation  time  In  samples. 

Hence,  under  this  assumption,  M  In  Eq.  (6)  Is  approximately  a  constant. 

Now.  make  the  above  reasonable  assumption,  that  M  Is  a  constant  over  the 
observation  time,  and  also  the  assumption  made  In  the  previous  section,  that 
the  time  samples  are  Independent  from  sample  to  sample.  Then,  the  probability 
density  function  for  the  N  samples  under  hypothesis  Hq  Is 

PqIv)’  PoU(l).  v(2).  ....  v(N)) 

N 

■  n 

n=l 

N 

-  ^  **(n)  M"’  v(n)  (7) 

The  exponent  In  Eq.  (7)  Is  expressible  In  terms  of  the  trace  function 


Tr{A)  = 


k 


E 


a 


11 


of  a  k  X  k  matrix,  A  *  Since  this  exponent  Is  a  scaler,  and  since 

Tr(A  B)  =  Tr(B  A),  where  A  and  B  are  matrices. 


(8) 


Hence,  In  terms  of  Eq.  (8),  Eq.  (7)  becomes 


-6- 


where 


-  Tr 


v(n)  v*(n)| 


-  N  Tr(H-\) 


(9) 


M  _ 

Mq  =  jr  v(n)  /(n)  =  V  V 

n=l 


(10) 


Is  the  well  known  MLE  of  the  unknown  covariance  matrix  under  hypothesis  Hq. 

A 

Here,  the  bar  operation  over  the  last  expression  for  Mq  denotes  the  time  average 

A 

over  n.  Note  that  the  Inverse  of  Mq  exists  only  If  K<  N. 

By  Eq.  (1),  for  the  signal -plus-noise  hypothesis  Hp 

E{y(n)lHi|  =  E{NQ(n)|  +  sjb  s(n)}=  b  s(n)  . 

Thus,  for  hypothesis  the  joint  probability  density  Is  found  in  a  similar 
fashion  to  be 


where,  by  Eq.  (1), 

-bl”)  “  ^  ’'l(")’  •••’  ’'k(")]^ 


(11) 


(12) 
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m 


i 


S  2  ■!!b(")V(") 

n=l 


-"b^b 


is  the  MLE  of  the  unknown  covariance  matrix  M  under  hypothesis 

The  parameter  vector  9  for  the  MLR  test  in  Eq.  (3)  consists  of  b,  the  com¬ 
plex  signal  envelope,  and  the  (k  +  2)(k+l)/2  functionally  independent  elements 
of  the  unknown  positive  definite  covariance  matrix  M.  It  is  more  convenient 
to  express  6  as  the  pair  ^  =  [b,  M]  so  that  the  sample  space  is 

=  {«}=  {[b,  M]|m  >o}  .  (14] 

where  M  >  0  denotes  the  fact  that  M  is  Hermitian  symmetric  and  positive  definite. 
With  this  definition  of  sample  space,  the  sets  w  and  A  -  w  associated  with  Hq 
and  H.J  in  Eq.  (2)  are,  respectively, 

w  =  {[0,  M]|m  >  o] 

A-  w  =  {[b,  M]|m>0,  b  oj  .  (is: 

By  Eqs.  (10),  (14),  and  (15),  the  MLEs  of  6  are 

©0  “  CO,  Mq]  and 

=  [b,  M^]  .  (is: 

with  respect  to  hypotheses  Hg  and  Hp  respectively,  holding  b  fixed.  A  substi¬ 
tution  of  Eq.  (16)  into  Eqs.  (3)  and  (4)  yields  the  MLR  test 


max  1; 


L(v)  = 


i  IlM  - 


mi  nil-  II N 
b  ”b 


>  K  then 
<  K  then  Hg  . 


This  test  is  evidently  equivalent  to 


misssmmm 


I 


je(v)  ■ 


>  C  then 
<  C  then  Hg  , 


(17) 


where  C  =  and  Mg  and  are  the  matrices  in  Eqs.  (10)  and  (13).  The  test 
in  Eq.  (17)  is  obtained  by  the  same  means  that  E.  J.  Kelly  [4]  used  recently  to 
find  adaptive  CFAR  tests  for  a  radar  signal  which  is  assumed  to  occur  in  only 
one  range  interval.  In  the  present  adaptive  CFAR  test,  the  radar  signal  is 
assumed  to  extend  over  many,  in  fact  N  »1,  range  intervals.  This  new  test  and 
the  methods  used  to  analyze  it  seem  to  differ  somewhat  from  the  Kelly  tests. 

The  test  in  Eq.  (17)  is  simplified  by  getting  rid  of  its  apparent  depen¬ 
dence  on  two  determinants  of  covariance  matrices.  This  is  to  be  accomplished 
by  the  well-known  Shur  identity  (e.g.,  see  [5,  pp.  45-47])  for  factoring 
determinants  of  partitioned  matrices  as  follows; 


A  B  II 

A  -  B  D"^  C 

0 

C  D  1 

C 

D 

s 

A  -  B  D“^  C 

D 

where  the  square  submatrix  D  is  assumed  to  be  nonsingular.  Thus,  if  N  >  k,  by 
Eqs.  (12)  and  (13)  and  the  Shur  identity. 


Hence , 


ll^bll  *(  |y  -  b  s|^  -  (y  -  b  s)  X*  x  (y  -  b  s)*  j||?^  || ,  (18a) 

A*  ■ 

where  it  is  convenient  to  let  M  *  x  x*.  Similarly,  by  Eqs.  (5)  and  (10), 

IKII  ^  (  |y|^  "  y  X  y*)  II  M  II  .  (18b) 

Substituting  Eqs.  (18a)  and  (18b)  into  the  ratio  test  of  Eq.  (17)  yields 

-  y  X  M  '  X  y 

J{v)<-~ - -  -  - — - - - ■■  -  >  C  for  H, 

( ly  -  b  s|^  -  (y  -  b  s)  x*  x  (y  -  b  s)*) 

“  '  "  <  C  for  Hq  .  (19) 

An  expansion  of  the  denominator  in  Eq.  (19)  produces 

D  =  ]yj^  -  b*  Ty*^  -  b  7^  +  b  b*  "[sl^ 

5F  -'.I  W  5F  ^.1  W 

-yx  M  £y+bsx  M  xy 

*  *  w  *  *  * 

+  byx^  M  £s-bbs)c  M  >(s 

=  C-  bB-b*B*  +  bb*A, 

where 


1  iZ 

w 

^-1 

M  ' 

- 9F 

A  = 

|s|  -  s 

X 

X  S 

* 

“T 

^-1 

M  ' 

- *• 

B  = 

s  y  -  s 

K 

X  y  ,  and 

C  = 

1  >> 

X 

''-1 

M  ' 

- w 

X  y  • 

* 

A  completion  of  the  square  in  the  above  quadratic  form  in  b  and  b  yields 

D  =  c  +  (Va  b  -  bYva j  I Va  b*  -  b^Va)  -  IbiYa 

=  C  +  \i7b  -  b^VaI^  -  >  0  . 


Evidently  D  is  minimized  when 


Va"  b  -  B*/  =  0  or  b  =  B^A 


so  that 


“n-in  '  -  I'ly*' 


V  ir  ^  » 

=  ly|  -yx  Mxy- 


s  y  -  s  X  M~  X 


Isl^  -  s  X*  X 


Substituting  Eq.  (20)  into  the  denominator  of  Eq.  (19)  gives 


2  ^  "I 


<  C  for  H, 


Clearly  this  test,  Eq.  (21),  is  equivalent  to  the  ratio  test. 


r(v)  = 


*  *■  -'-l  2 

sy-sx  M 


*  i 

s  X  M  X  s 


S-'  xy*) 


>  for  H, 


<  '"o  *’»'•  "o  •  < 


where  r^  =  1  -  1/C. 
0 


It  is  now  shown  that  the  range  of  ratio  r(^)  in  Eq.  (22)  is  restricted 


to  the  interval  between  zero  and  one,  i.e.  0  <  r(v)  <  1.  In  order  to  show  this 


it  is  convenient  to  let 


e  -  X 


(23) 


/2  A. 

where  M  '  is  the  inverse  of  the  square  root  of  M.  For  any  given  _x  or  M  = 

2^  >(  ,  the  square  root  may  be  defined  by  M  '  =  U  A  '  U  ,  where  A  '  is  diagonal 

A 

with  elements  equal  to  the  square  roots  of  the  eigenvalues  of  M  and  U  a  unitary 
matrix.  With  these  definitions,  one  has  the  identity 


XX*  frV2 


=  H" 


1/2  ii  ft-1/2. 


Ii,  . 


(24) 


where  I|^  denotes  the  k  dimensional  identity  matrix,  and  the  upper  bar  operation 
again  denotes  the  average  over  the  N  samples. 

The  key  identity  needed  to  find  the  range  of  r(v)  is  the  identity 


s  y  -  s 


7y^  *  [s  .  (s  /)  e]  [y  -  (y7^)  e]' 


(25a) 


To  prove  Eq.  (25a),  use  Eq.  (24)  as  follows; 


If  one  lets  y  =  s  and  s  »  y  in  Eq.  (25a),  the  following  corollary  identities 
obtain: 


and 


(2! 


That  the  range  of  r{^)  is  0  <  r{^)  <  1  follows  from  Eqs.  (22),  (25a),  (25b),  and 
an  evident  application  of  Schwarz  inequality  to  the  right  side  of  Eq.  (25a),  as 
follows: 


[s  -  (s^)  e][y  -  (7i*)  e] 

<  Is  -  (7^)  el^  ly  -  {y^)  ej^  . 

IV.  FALSE  ALARM  PROBABILITY  FOR  CFAR  SIC  TEST 

In  order  to  evaluate  performance  of  the  CFAR  SLC  test  derived  in  the  pre¬ 
ceding  section,  it  is  important  to  estimate  the  probability  of  a  false  alarm  as 
a  function  of  detection  threshold  r^,  i.e.,  one  needs 

PFA  =  Prob  |r(^)  >  t^IHq  |. 

Hence,  the  probability  density  function  (PDF)  of  the  test  function  r(v)  must  be 
obtained  in  order  to  compute  the  PFA. 

To  derive  the  PDF  of  r(\^),  it  is  convenient  for  preciseness  to  perform  the 
averaging  operations  which  appear  in  r(v)  as  vector  operations.  This  is  accom¬ 
plished  by  changing  the  functions  s  =  s(n),  y  =  y(n),  *  x{n)t  and  e  =  e(n) 

for  1  <  n  <  N  into  vector  and  matrix  notation,  as  follows: 
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^  ... 


) 


yi 


3 


H.. 


•S4? 


s  =  [s(l),  s(2).  ...,  s(N)]  ,  a  row  n -vector 

Z  *  [y(l).  y(2),  ....  y(N)]  ,  a  row  n -vector 

X  =  [x(l)  :  x(2)  ;  ***  :  x(N)].  a  k  x  N  matrix 

J  =  [e(l)  :  e(2)  :  ***  :  e(N)|,  a  k  x  N  matrix  . 


In  terms  of  this  new  notation. 


X  X*  =  N  M  and  J  *  (x  X*)"^/^  X  . 


Note  in  this  new  notation  that  the  identity  in  Eq.  (24)  corresponds  to 


JJ*-(x  X*)-’/2x  X*  (x  X*)-’'2  .  I 


(27a) 


and  that  r  in  Eq.  (22)  is  expressible  as 


r(v)  = 


Is/  -  s  x*(x  x*)-' _ 

(s  s*  -  s  X*  (X  x’)-'  X  £  )(££  -  JL  X*  (x  **)■’  »  /) 


*  *  *  2 
s_^  -  s_0  0  ^ 


I  fr  *  Wv  /  ♦  h  *v 

ill  “1'^  ‘^ijiy.^  -  y.  J  J  y.  )  . 


(27b) 


Next,  note  the  identity 


j)  (jr  -  U  0  j)‘ 

-  Q,(£)  *  ( 

where  Q-j  and  Q2  are  quadratic  forms  in  the  components  of  (see  Appendix  A). 
This  identity  is  the  same  as  the  identity  in  Eq.  (25b),  but  in  the  notation  of 
Eq.  (26).  It  follows,  by  expanding  the  last  term  of  Eq.  (28)  and  using  Eq. 
(27a). 


From  (26)  and  (27),  the  rows  of  J  are  linearly  independent— in  fact,  ortho- 
normal  in  an  inner  product  sense.  Also,  by  Eq.  (26),  the  Hermitian  symmetric 


matrix  J  J  can  be  expressed  as 


* 

J  J  = 


e*(l)' 

e*(2) 


Vd)  j' 

e*(2)  J 


e*(N) 


e*(N)  J 


so  that  the  i-th  row  of  J  J  is  ^  (1)  J.  Hence,  every  row  of  J  J  is  a  linear 
combination  of  the  k  rows  of  J.  Since  the  k  rows  of  J  are  linearly  Independent, 

•k 

the  rank  of  J  J  is  k. 

For  any  given  X  or  J,  a  unitary  transformation  can  be  found  which  removes 
the  direct  functional  dependence  of  the  ratio  r  in  Eq.  (27b)  upon  either  of 
these  matrices.  The  existence  of  such  a  unitary  transformation  is  guaranteed 
by  the  generalization  of  Cochran's  theorem  given  in  Appendix  A. 

Consider  any  N  x  N  unitary  matrix  U  of  form 


“k+1 ,1 »  • • • *  “k+1 ,N 
U  = 


u»  1 ,  . , . ,  u. 


and  the  unitary  transformation 


(29a) 


(29b) 


This  transformation  first  sends  the  left  side  of  Eq.  (28)  into 


*lUU  £  =££  «  2^  hj  ' 

1=1 

It  next  transforms  the  first  quadratic  form  on  the  right  side  Into 

<ii(y)  '(ji  01/ 0  =  h . 9k][9, . gj* 

- 1  • 

1=1 

Thus,  by  the  same  transformation,  the  second  quadratic  form  Is  transformed 
Into 

^2(y)  "2  hjp  ■  53  hjp '  53  ho  ^  * 

1=1  j=l  j=k+l 

so  that  the  rank  of  Q2(j!l)  N  -  k,  I.e.,  It  has  (N  -  k)  non-zero  eigenvalues  and  k 
eigenvalues  equal  to  0.  Since  the  requirements  of  Cochran's  theorem  are  satis¬ 
fied,  a  unitary  matrix  of  the  form  of  Eq.  (29)  exists  In  actuality,  which  will 
transform  the  second  quadratic  form  Into  a  sum  of  squares  In  such  a  manner  that 
Ql  and  Q2  have  no  variable  In  common. 

Now,  apply  this  unitary  transformation  to  the  ratio  r  in  Eq.  (27b)  and,  as 
well,  transform  the  signal  £  by  letting 


Then,  the  ratio  In  Eq.  (27b)  evidently  becomes 


^tuu  t  |tj|-Hauu  a  -  ^ 

I  "  IP 

L  •/ 

*7"TI — - T  ■ 

(E  l•.ll  E  \;f) 

\j=k+l  /\)*k+l  / 

For  any  particular  J,  as  defined  In  Eq.  (26),  or  set|x(n)|n  *  1,  2,  ..., 
njof  AUX  data,  the  unitary  matrix  U  In  Eq.  (29a)  and  finally  the  vector  £ 
from  ^  In  Eq.  (29b)  can  be  determined.  Thus,  for  the  no1se-on1y  hypothesis  Hq 
and  for  a  given  J  or  X,  the  resulting  components  of  £  are  mutually  statistically 
Independent.  This  follows  from  the  assumption  made  In  Section  II,  that  the 
samples  of  y(n)  were  zero-mean  Gaussian  and  mutually  Independent,  and  from  the 


following  argument: 


e{i*  ajx.  «o}  °  “ 


=  U  E 


=  2  OrM, 


i  ai)*!*.  "o) 
{a*  £|Ho}  U*  •  u( 


U*  =  u(2<r2 


where  Is  the  noise  power  per  sampled,  quadrature  demodulated,  component  of 
y(n)  for  1  <  n  <  N. 

Now,  In  Eq.  (30),  let 


w..  ‘  tJ 


A/ V 
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m 


a 


r- 


■ij 


Ui 


for  k+1  <  j  <  N,  then  Eq.  (30)  simplifies  further  to 


E  "jfj* 

j-k+1 


j»k+l 


Another  Identity  similar  to  Eq.  (28)  can  now  be  used  to  find  another  uni¬ 


tary  transformation  to  eliminate  the  functional  dependence  of  r  on  the  w.  In 

J 


Eq.  (31).  This  Identity  Is 


N  N  2  H  /  H  \* 

E  kip  •  E  “i  9i*  *  E  -( S  "t  Si*)" 

1=k+l  1»k+l  1*k+l  \l=k+l  / 


*  12 


“  ^3K+1*  •••*%)■*■  Q4K+I 


The  two  quadratic  forms  and  of  this  expression  can  be  treated  In  the  same 
manner  as  In  Eq.  (28)  to  show  that  has  rank  1. 

Next,  consider  any  (N  -  k)  x  (N  -  k)  unitary  matrix  of  form 


.  .  .  ,  W|y 


^k+2,  k+1,  ^k+2,  k+2,  **k+2,  N 


(33a) 


^N,  k+1  *  ‘‘n,  k+2  . ‘^N,  N 


and  the  unitary  transformation 


Zi,  .  1  ,  . . .  ,  Z| 


1]  •  kk+i 


* 

D  . 


(33b) 


This  transforms  the  left  side  of  Eq.  (32)  Into 
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E  hr 

j=k+1 

and  Qj  into  Hence,  is  transformed  into 

N 

E  r.r 

j-k+2 

and  has  rank  N  -  k  -  1.  Again,  the  requirements  of  Cochran's  theorem  are 
satisfied  and  a  unitary  matrix  of  form  0  exists  which  will  transform  both 
and  into  a  sum  of  squares  with  no  variables  in  comnon. 

If  one  applies  the  unitary  transformation  of  Eq.  (33b)  to  ratio  r  in 
Eq.  (31),  one  obtains 


r(v)  •  -  (34a) 

E  h/ 

J-k+1 

as  the  test  function  on  the  condition  that  X  or  x(n)  for  (n  ■  1,  2,  ...,  k) 
is  given.  By  the  same  procedure  as  before,  the  z.'s  are  mutually  independent. 

J 

But,  also  by  Eq.  (1),  under  hypothesis  Hq,  the  quadrature  components  of  y(n) 
and  each  component  of  x(n)  are  independent.  Furthermore,  a  unitary  transforma¬ 
tion  on  a  complex  unitary  space  is  equivalent  to  an  orthogonal  transformation 
of  twice  the  dimension.  Thus,  the  numerator  of  Eq.  (34a)  consists  of  two  squares 
of  independent  real  variables  and  the  denominator  consists  of  2  (n  -  k)  squares 
of  mutually  independent  real  variables.  Hence,  the  ratio  in  Eq.  (33a)  has  the 


r'(v)  «  — ^ 


R  2  +  T  2 

•^k+l  *k+1 


E  («/  ^ ./) ' 


(34b) 


■riyfi 


Zj  »  Rj  +  1  Ij  for  (j  »  k+],  ....  n) 
and 

e[r/|x,  H„].e[i/|x.  Ho]-<r2. 

The  form  of  the  conditional  random  variable  r(v^)  in  Eq.  (33a)  or  (33b) 
suggests  that  the  conditional  probability  density  of  r(^)  is  a  generalized 
Student's  t  or  Beta  density  (see  Ref.  2.  pp.  237-243).  Thus,  by  Ref.  2 
(18.4.2), 

Prob  |r  <  r(v)  <  r  +  d  »*jHQ,  x|  =  P(r|HQ,  x)  d  r 


where,  for  the  present  case,  m  *  2  and  n  =  2  (N  -  k)  -  2.  Hence, 

'’(••K- X)  ■  TtMt 

-  (N  -  k  -  1)  (1  -  r)^"*^"^  . 

Since  the  right  side  of  Eq.  (35)  is  functionally  independent  of  X  or 
2(^(n)  for  (n  =  1,  2,  ...,  k),  the  absolute  probability  density  of  rjy)  under 
hypothesis  Hq  must  be  the  same,  i.e., 

p(r|HQ)  =  (N  -  k  -  1)  (1  -  r)'^"’^"^  . 

Integrating  P(r|HQ)  from  the  detection  threshold  r^  to  1  yields 


PFA  »  Prob  |r(v)  >  i^qIhqI 

•  j  (N  -  k  -  1)  (1  -  dr 


=  (N  -  k  -  1)  I- 


=  ('  -  '•o) 


N-k-1 


for  the  actual  probability  of  a  false  alarm  for  the  CFAR  sidelobe  canceller  cri¬ 
terion  found  in  Eq.  (22). 


V.  CONCLUDING  REMARKS 

A  general  CFAR  signal  detection  test  has  been  found  for  a  k  -  AUX  sidelobe 
canceller  in  Eq.  (22).  For  k  -  0,  this  test  reduces  to  the  classical  CFAR  power 
ratio  test  (see  Ref.  3,  p.  326,  for  the  first  such  test— the  normalized  periodo- 
gram  test). 

In  Section  IV,  the  formula  for  the  probability  of  a  false  alarm  is  found 
and  given  in  Eq.  (37).  When  k  -  0,  this  result  reduces  to  the  false  alarm  rate 
given  for  the  usual  CFAR  test.  The  methods  used  to  find  the  PFA  for  the  test 
in  Eq.  (22)  are  similar  to  the  techniques  Fisher  used  to  find  the  Fisher  z  - 
distrubution,  e.g.,  see  Ref.  2,  sections  11.11,  29.1,  and  29.2. 

For  a  complete  performance  study  of  the  CFAR  SLC  test  developed  in  this 
report,  it  would  be  desirable  to  find  a  formula,  possibly  approximate,  for 
the  probability  of  detection.  It  is  expected  that  it  will  involve  the  tech¬ 
niques  developed  and  used  in  Refs.  6  and  4. 


REFERENCES 


1.  Brennan,  L.  E.,  and  I.  S.  Reed,  A  Constant  False  Alaam  Rate  Sidelobe  Can¬ 

celler  for  Radar,  Third  Quarterly  Progress  Report,  submitted  to  the 
Naval  Air  Systems  Command  on  Contract  NOOOl 9-84-C-0074  by  Adaptive  Sen¬ 
sors,  Inc.,  February  1985. 

2.  Cramer,  H.,  Mathematical  Methods  of  Statistics,  Princeton  University  Press, 

1963. 

3.  Kelly,  E.  J.,  I.  S.  Reed,  and  W.  L.  Root,  "The  Detection  of  Radar  Echoes 

in  Noise,  Part  I,"  Society  for  Industrial  and  Applied  Mathematics,  8(2), 
June  1960,  pp.  309-341. 

4.  Kelly,  E.  J.,  An  Adaptive  Detection  Algorithm,  a  forthcoming  Technical 

Report  from  MIT  Lincoln  Laboratory,  Lexington,  MA. 

5.  Gantmacher,  F.  R.,  Matrix  Theory,  Vol.  1,  Chelsen  Publishing  Co.,  New  York 

1977. 

6.  Reed,  I.  S.,  J.  0.  Mallet,  and  L.  E.  Brennan,  "Rapid  Convergence  Rate  in 

Adaptive  Arrays,"  IEEE  Transactions  on  Aerospace  and  Electronic  Systems, 
AES-10,  November  1974,  pp.  853-863. 


Appendix  A 

UNITARY  TRANSFORMATION 


Let  i  \ 
with  elements  £  =  »  x 
x^  of  form 


be  the  N  components  of  an  N  dimensional  unitary  space 
2*  •••>  xJ.  Consider  quadratic  forms  in  x-i*  X2f  ... 


N 

oh-  *2 . "n)*  2  "ikN’'/- 

i,k=l 

where  h|^^-  =  and  where  the  asterisk  denotes  complex  conjugate.  If  Q(xi  *  ... 
x^j)  >  0,  Q  is  a  non-negative  quadratic  form.  If  this  relation  equals  zero  only 
when  all  the  x^  equal  zero,  Q  is  said  to  be  positive  definite.  The  above  matrix 
H  ■  with  the  property  H  =  H*  is  called  Hemitian  symmetric.  If  the  rank 
r  of  H  is  less  than  N,  H  and  its  corresponding  quadratic  form  Q  are  called  posi¬ 
tive  semi -definite.  If  rank  of  H  is  r,  the  rank  of  its  quadratic  form  Q  is 
also  r. 

The  generalization  of  Cochran's  theorem  [2,  Sec.  11.11]  needed  to  find 
the  PDF  of  r(^)  in  Eq.  (27b)  of  Section  IV  is  as  follows: 


Theorem:  Let  an  identity  of  form 


N 

l^il  ”  ^1  ^  ^  ^t 

i=l 

be  given,  where  Q^.  for  i  *1,2,  ...,  t  is  a  non-negative  quadratic  form  in 
X, ,  x,,  . . . ,  Xu,  of ,  at  most,  rank  r. .  Then,  if 
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there  exists  a  unitary  matrix  U  and  its  corresponding  unitary  transformation 
X  =  U*,  which  changes  each  Q^.  into  a  sum  of  the  squares  of  the  absolute 

values  of  the  y^. ,  according  to  the  relations 

■•l  N 

'*1*2  |^l|^  ■  ‘>2  “  S  . ‘’k  "  2  |^l|^  ’ 

1  r,^.,  N-r^+1 

in  such  a  manner  that  no  two  Q.  have  a  variable  y.  in  common  and  where  x  = 
j^Xp  ^2*  ....  X|^ j  and  y  =  ^y-j ,  y2.  ...»  yj^j  are  row  vectors  and  "**'  denotes 
conjugate  transpose. 

This  proof  of  this  theorem  is  precisely  the  same  as  the  proof  given  by 
Cramer  [2],  except  that  wherever  a  square  of  a  quantity  appears  in  his  proof 
it  is  here  a  square  of  an  absolute  value.  Also,  replace  the  words  "orthogonal 
transformation"  in  the  Cramer  proof  by  "unitary  transformation."  Every  other 
step  and  statement  in  the  proof  remain  the  same. 


